Interacting trapped bosons yield fragmented condensate states in low dimensions 
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We investigate the level population statistics and degree of coherence encoded in the single-particle 
density matrix of harmonically trapped low-dimensional [quasi-one-dimensional (quasi- ID) or quasi- 
two-dimensional (quasi-2D)] Bose gases with repulsive contact interactions. Using a variational 
analysis, we derive fragmentation of the condensate in the weakly confining directions into two 
(quasi-lD) respectively three (quasi-2D) mutually incoherent macroscopic pieces, upon increasing 
a dimensionless interaction measure beyond a critical value. Fragmented condensate many-body 
states in low-dimensional systems therefore occur well before the thermodynamic limit of infinite 
extension is reached, in which phase fluctuations of the matter wave field create an infinite number 
of nonmacroscopic fragments. 

PACS numbers: 03.75.Gg 



I. INTRODUCTION 

Bose-Einstcin condensation, the macroscopic occupa- 
tion of one field operator mode (= single-particle orbital), 
ceases to exist in the thermodynamic limit of infinite ex- 
tension in spatial dimension less or equal to two [ll, ^, 
as a direct consequence of the Bogoliubov inequality 3| . 
This so-called Hohenberg-Mermin- Wagner theorem has 
been argued also to hold for systems with a nonvanishing 
cross-section (quasi-lD) or a finite thickness (quasi-2D) 
[1]. On the other hand, for a finite extension respec- 
tively radius of curvature along the cylinder axis or in 
the plane, condensation can persist Q. Low-dimensional 
Bose gases are nowadays routinely created by varying 
cloud aspect ratios in strongly anisotropic trapping ge- 
ometries of magneto-optical or purely optical origin [6|-l8| . 
While the very existence of condensates does not depend 
on interaction in the thermodynamic limit, the increasing 
dependence of their detailed dynamical and static proper- 
ties on the interaction coupling when lowering the dimen- 
sionality has been demonstrated experimentally 0, [l^ as 
well as investigated theoretically [ill, Ell • 

The search for fragmented condensate states of Bose 
gases, that is the macroscopic occupation of more than 
one field operator mode with no mutual coherence be- 
tween them, has a long history, see, e.g., [isl - fisj . The 
obvious general difficulty in answering the rather deli- 
cate question whether a given system is coherent (a con- 
densate) or fragmented consists in solving an interacting 
many-body problem, from which the single-particle den- 
sity matrix follows. The latter, by definition [l^, gives 
the degree of fragmentation versus that of any remaining 
coherence. 

In the following, we shall show that fragmentation of 
condensates takes place well before the thermodynamic 
limit is reached in lower-dimensional Bose gases. We thus 
argue that interaction-induced fragmented phases, with 
only a few macroscopically occupied modes (termed frag- 



mented condensates [l6|). exist before fragmentation into 
infinitely many incoherent pieces of nonmacroscopic size 
sets in, occurring when in the thermodynamic limit phase 
fluctuations destroy the condensate in the quasi-lD or 
quasi-2D low-dimensional regimes [l|, 0] . 

We concentrate on the point where fragmentation sets 
in first, into two (quasi-lD) respectively three (quasi- 
2D trapping geometry) modes with no coherence be- 
tween them, which may be described by the solution of a 
many-body problem essentially exactly solvable for any 
strength of interaction. Thus we demonstrate how the 
formalism laid down in [l7| leads to fragmentation of a 
scalar Bose gas in the weakly confining directions. It is 
stressed that for a trapped system positive interaction 
couplings lead to fragmentation above a critical value, 
counter to the expectation from the nontrapped thermo- 
dynamic continuum limit, where the physics is homoge- 
neous and the Fock exchange term in the total energy is 
preventing fragmentation [14| . 



II. QUASI-ID TRAPPING GEOMETRY 

We begin with the simple case of one excitation direc- 
tion along the weakly confining (z) axis of a Bose gas in 
the quasi- ID limit. Assuming that only two (the ener- 
getically lowest) longitudinal modes are significantly oc- 
cupied [that is, with macroscopic occupation numbers of 
0{N)], we restrict the expansion of the field operator to 
these two modes, ^(r) = X]i=o i the modes are 

normalized to unity, / <Pr\^i{r)\^ = 1. 

The validity of the two-mode approximation on the 
mean-field level and in double-well traps has been dis- 
cussed in [3 [111 ■ the following, we go beyond mean- 
field (solve for the full quantum solution) and consider 
a single trap, determining the interaction-dependent 
orbital parameters self-consistently by looking for the 
ground state variationally (see below). The two-mode 



2 



approximation can then be expected to represent accu- 
rately the interacting many-body physics in a quasi-lD 
trapping setup for much larger values of the coupling 
constant than one would anticipate from a mean-field 
analysis with fixed orbitals [Toj . 

Consider the interacting two-mode Hamiltonian 
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+ —a{a{aoaa + —a^al^aiai + —a^aia'^aa, (f) 



with the interaction coefficients Ai = Vbooo,^2 = 
V^iiii,^3 = Viioo,^5 = ^0011, and A^ = 1/0101 + Vww + 
^1001 + ^110 1 where the interaction matrix elements are 
given by V^jki = g ^ (Pr-^*{r)-^*{r)-^kir)-^i{r). The 
contact interaction is assumed to be repulsive, g = 
AttUs > 0, with fls the s-wavc scattering length (we put 
h ~ m — 1^ where m is the boson mass); the single- 
particle energies ei > eg- The pair-exchange coefficients 
A3 are real or complex numbers depending on the choice 
of modes {Ai, A2, A4 are always real by definition), and 
in general do not vanish in a spatially confined gas. The 
spatial basis dependence of the Hamiltonian stemming 
from the position dependence of the modes' phase is con- 
tained in the values of the A3, and is reflected in the 
correlation functions characterizing the response of the 
system to external perturbations. Global phase shifts, 
^k{r) — ^k{r)e'^^^ with 9k constants, as is well known, 
leave the correlation functions invariant. 

We perform an expansion of the two-mode many-body 
wavefunction in a Fock basis. 
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so that ipi is the probability amplitude_ 
tides residing in the excited state 
total density then is given by {'^\p{r)\^) 



Nonvanishing pair coherence 
expectation value J^/it/i^i^i^ 
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as defined by 
2 x^*! OQaJSiai + a|a|aoao|*) = 

where the pair-exchange coefflcient 



di = ^{l + 2){l + 1){N - I - l)(iV - I), is enforced by 
energy minimization, and yields definite a value of A3 
and coherent pair oscillations between the two modes, 
with A3 equal to their frequency. These coherent pair 
oscillations are thus in analogy to Josephson oscillations, 
which are due to an off-diagonal matrix element of the 
single-particle density matrix (\E'| aJSi \'^) oc N, analo- 
gously represented by single-particle exchange terms of 
the form — ^ajai + h.c. in the Hamiltonian, where il is 
the Josephson frequency. 

Minimizing the energy with respect to the Fock state 
amplitudes 0/, the linear system to be solved reads 



A, 



{dl'lpl+2 + dl-2'ipl-2) + Cllpl 



(3) 



where E = is the total energy. The diagonal 

coefficient reads q eo(iV - I) + + ^Ai{N - 1){N - 
l-l) + \A2l{l - 1) + \Ai{N - 1)1. The above discrete 
equation connects the values of V'z on Z- "sites" differing 
by two, and leads to alternating signs within the even and 
odd / sectors of the -0; when A3 > 0, sign(?/;i^i+2) = ^ 1- 
As a result, sign{ipiipi+i) = ±(—1)', the ± sign reflecting 
the initial condition on the ipi [20|. This, then, leads to 
the destruction of first-order coherence defined by gi = 

\ (4-1 ajai + ajao |*) = Re [E^Ii '/'f-i^i - « + 1)] . 

provided A^^ A2-V 2 A3 - Aa>^ 

We take as the (real) modes the ground and first ex- 
cited state oscillator modes along the cylinder axis 
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where Rz is a variational parameter for the orbital de- 
localization length away from the harmonic oscillator 

— 1/2 

ground state in which R^ = Iz ^ ^z, wj^ are the 

harmonic trapping frequencies in the z and radial direc- 
tions, respectively, where Iz = oJz ^^"^ ■, l± = ^±^^'^ are the 
corresponding harmonic oscillator lengths. The trans- 
verse degrees of freedom are assumed to be frozen, so 
that the transverse part of the wave function is the har- 
monic oscillator ground state. The above choice for the 
modes leads to A3 real and positive (note that generally 
A4 = 4A3 for real modes and contact interaction). 

The single-particle energies of the motion along z, oc- 
curring in the Hamiltonian (jlJ [the single-particle en- 
ergy in the transverse direction is dynamically irrele- 
vant], are eo = jR^uJl + and ei = Seo, while the 
interaction coefficients for these orbitals are Ai = Uq = 
g/{Rzll{2TTf'^), A,/Ai = {f,|,i,2}. Analytical ex- 
pressions for the ground-state distribution can be ob- 
tained when the absolute value of the Fock state am- 
plitudes \ is considered as a continuum variable [2ll |. 
In the continuum limit we have a Gaussian distribution 
for the absolute value of 
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where 



the 
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A1+A2+2A3-A4 
state centered at I 



effective oscillator length a^^^ = 
and the shift from a fully fragmented 
reads 6 



N_ 
2 



N{Ai-A2)/2+eo-ei 



[l7|. The degree of fragmentation is in the two- mode 
model defined hy ^ ^ f — |Ao — )^i\/N, where Ao,i are 
the two eigenvalues of the single-particle density matrix 
resulting in {Ni = (n.^) = (aj'ai)) 



^l-±i^NoN,-\{dldrW). (6) 
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Because in the present case first-order coherence (aJSi) 
essentially vanishes [lO], the fragmentation measure j 
is simply determined by the mean occupation numbers 
of the two states. In the continuum limit, the degree 
of fragmentation is obtained to be = 1 — -^l^. Two- 
mode fragmentation is therefore in this limit consistently 
obtained when Oosc is real as well as © < N /2 holds, and 
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is becoming maximal when © = 

The total energy is in the continuum limit ([5]) for the 
two-mode model ([T]) given by 
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The relative values of the interaction coefficients Ai lead 



to ai 



|iV, which results in = y + - 



^X^. The ratio X = eo/(iVAi) = ^ + A^), is a 
measure of the ratio of single particle and interaction en- 
ergies, where the dimensionless variational parameter is 
defined by A = Rjl^. The shift reads © = f (1 - 16X), 
and the quantity 



Gi 



(8) 



is a dimensionless measure of interaction strength. We 
used in the final expression on the right-hand side 
that the quasi-lD interaction strength reads giD ~ 
f//(27r^^) = 2as/l\ (valid when Og <SC i.e. away from 
geometric scattering resonances [lH), with dimension of 
inverse length. Note that Gi linearly increases with in- 
creasing Iz, indicating that single condensates cease to 
exist for Iz going to infinity for any particle number 
and at any finite interaction strength, which is in accor- 
dance with what we would expect from previous studies 

The total energy obtained from solving the discrete 
many-body equations is to be minimized to obtain A = 
A(Gi). The analysis reveals that the numerically deter- 
mined energy as a function of A, for a given Gi, has 
only one, nonfragmented, minimum for small enough Gi 
below a critical value, while above the latter, a second 
local minimum appears. The transition to a fragmented 
state, corresponding to the second minimum at a smaller 
(scaled) delocalization length A, is discontinous, and the 
degree of fragmentation jumps to a finite value, cf. Fig.[T] 
Therefore, the phase transition between single conden- 
sate and fragmented condensate states obtained within 
the two-mode model ^ is of first order. Increasing the 
dimensionless interaction measure further, the nonfrag- 
mented local minimum gradually disappears (asymptot- 
ically changing into a turning point), and only the frag- 
mented minimum of the total two-mode energy remains 
(top curve in Fig.[T]). 

We obtain from the continuum limit Eq. ^ for the 
energy per particle in terms of Gi and A the following 



FIG. 1. (Color online) Evolution (not to scale) of the energy 
landscape as a function of the variational parameter A, with 
increasing dimensionless interaction measure, from bottom to 
top [the energy minima have additionally been shifted rela- 
tive to each other for clarity] . Energy curves are from the full 
numerical solution of The continuum limit is indicated 
by the dashed curve for the energy curve close to critical- 
ity (second from the bottom), showing that the continuum 
limit energy, tending to negative infinity for large A, does 
not describe the nonfragmented minima of the total energy, 
cf. Eq. Beyond a critical interaction measure, a new mini- 
mum corresponding to fragmentation appears (empty square) , 
with a discontinous jump from the nonfragmented minimum 
(circle) to a smaller A (third curve from the bottom). In- 
creasing the interaction measure further, the nonfragmented 
minimum eventually completely disappears. 



polynomial equation, using that GiUJzl Ai) ~ A, 



E 



7/1 



1 

6Ch 



A3 



AM + ^ 



Gi 
3A' 

(9) 



For large A, we may approximate the above equation 
as ~ -I- The extrema equation is 

thus quadratic in A^/Gi, so that A is scaling with G^^^ 
in this limit; the minimum solution describing a frag- 
mented state is A^ = ^Gi (the other solution does not 
describe a minimum, cf. Fig.[l] dashed line). Note that 
the large Gi scaling A cx {Naslz/l^y^^ is identical to 
the single-condensate Thomas-Fermi scaling Q [with a 
different prefactor]; the quantity X = eo/{NAi) asymp- 
totes to 1/10 in this limit. 

Evaluating in the large Gi limit the difference be- 
tween the fragmented energy minimum of ([5]) at A = 
(2Gi/5)i/3 and the energy obtained from a single con- 
densate of all particles residing in the energetically lower 



mode, E, 
A 



gY^, we get H 
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which has a minimum at 



AE 



E, 
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0.02 X Gi^^ 



(10) 
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The energy difference between single condensate and 
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fragmented condensate ground states is therefore very 
small around the transition point, increasing with the 
particle number like N^^^. 

We obtain in the continuum limit for the degree of 
fragmentation 



S{Gi) = 1 
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where the last inequality represents the large Gi — large 
A limit. The maximally achievable two- mode degree of 
fragmentation from this variational ansatz for the modes 
is thus 80 %. We conclude that fragmentation of the con- 
densate many-body state occurs for sufficiently strong in- 
teraction coupling at a given finite extent Rz of the cloud, 
that is when 



Gi > (Gi)c^9.8 



(12) 



where the critical value (Gi)c for a fragmented minimum 
of the total energy to first appear, cf. Fig[Tl has been de- 
termined from the full numerical solution of ^ . We dis- 
play the continuum result for ^ = S'(Gi) obtained from 
the solution of Eq. in Fig. [21 The system is for small 
interactions and densities (Gi <C 1) a single condensate 
and A is close to unity; the slow increase of the optimal A 
with Gi additionally demonstrates that an ansatz incor- 
porating only the two lowest harmonic oscillator states 
is a reasonably accurate approximation. For large shift 
(large A), the continuum description fails, and the full 
discrete coupled system of equations for ipi in Eq. ^ has 
to be solved to find the optimal value of A. The result 
from the (exact) two-mode many-body wavefunction cor- 
responding to the numerical solution of ^ is displayed 
together with the continuum limit in Fig. [21 where very 
good agreement is visible. 



III. QUASI-2D TRAPPING GEOMETRY 

Consider now a cylindrically symmetric trap, with an 
infinite continuum of excitation directions in real space, 
that is, in all radial directions. A crucial difference to 
the quasi-lD case is the topology, resulting in an an- 
gular quantum number m^. Furthermore, the lowest 
(7710 = ±1) azimuthal excitation modes of the cylin- 
drically symmetric harmonic oscillator are lower in en- 
ergy than the radial excitations. The harmonic oscillator 
ground state in radial as well as transverse directions and 
the two lowest excited states along the azimuthal direc- 
tion are given by 
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and arc thus the three lowest-lying single-particle states. 
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FIG. 2. (Color online) Optimal value of the variational pa- 
rameter A, and degree of fragmentation 5 (inset), both as a 
function of the dimensionless measure of interaction strength 
Gi, in the quasi-lD trapping geometry. Numerical solutions 
of the discrete set of many-body equations ^ are depicted 
by blue circles (nonfragmented energy minima), and by red 
squares, respectively, the latter describing the energy minima 
at which fragmentation occurs, cf. Fig.[T] The solution in the 
continuum limit of Eq. ((Ojl is shown for the fragmented state 
by the black solid line. The number of particles is A'^ = 500. 



A. Two-mode approximation 

We first investigate whether fragmentation occurs be- 
tween the ground state and a (coherent) superposition of 
the left- and right-circulating azimuthal waves ip^ , ip- . 
We shall consider below the full three-mode problem of 
taking all three states in (jl3p into account, and will show 
that in the continuum limit the corresponding predic- 
tions become essentially identical to those of the simpler 
two-mode model, which we discuss here to most clearly 
elucidate the primary differences to the quasi- ID trap- 
ping geometry. 

We consider an equal- weight coherent superposition for 
the excited state wavefunction (such that it has angular 
momentum zero like the ground state). 



^^/^R\l, 



■ r cos (p exp 
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and evaluate the coefficients parametrizing the two-mode 
Hamiltonian We obtain Ai = g/{{2T:f/'^R\lz) = 
Vo = ^ and A/Ai = {l,f,i,2} with A = R^/lr, 
the ratios of the interaction coefficients are therefore 
identical to the quasi-lD case. We furthermore have 
i?^ and £1 = 2eo. The result for the 
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FIG. 3. (Color online) Energy minima of the two-mode ansatz 
for the quasi-2D trapping geometry in terms of the variational 
parameter A, and degree of fragmentation 5 (inset), both as 
a function of the interaction measure G2. Conventions are 
the same as those in the quasi- ID trapping setup of Fig. [2] 
Note that the decrease of the size A for fragmented states 
relative to that of single condensate states is significantly less 
than in the quasi-lD case. The number of particles here is 
A'^ = 20000; a higher particle number is necessary to achieve 
agreement of numerics and analytics comparable to that of 
the quasi-lD case. 



at the fragmented energy minimum. The large G2 scaling 
A cx (Nas/lzY^'^ is, like in the quasi-lD setup, identical 
to the single-condensate Thomas- Fermi scaling @; the 
quantity X = eo/{NAi) asymptotes here to 1/3. In the 
large G2 limit the energy difference between the frag- 
mented energy minimum of (|16p and the energy mini- 
mum of Egs = -I- A^) -I- G2/{2K^) [again assum- 
ing all particles to be in the energetically lower mode, 
i.e., that the many-body state is a simple condensate or 
equivalently a Z = Fock state], at A = G^^^ , is given by 

W^^-^-lir^^ 0.002 X G\'\ The energy differ- 
ence between condensate and fragmented ground states 
thus increases more slowly than in the quasi- ID case with 
the dimensionless interaction measure, cf. Eq. (|1C 
The two-mode fragmentation measure gives 
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The maximally achievable two-mode degree of fragmen- 
tation from the ansatz (|14p is thus 44%. Fragmentation 
takes place when the dimensionless coupling (cf. Fig.|3l 
where we compare numerical and continuum limit ana- 
lytical results for the fragmentation) 



G2 > (G2). 



80. 



(19) 



It involves g2D = g/W'2'n'lz), the quasi-2D coupling con- 
stant when fls, i-C, away from geometric scattering 
resonances [23. Observe that the dimensionless inter- 
action strength G2 here is independent of the harmonic 
oscillator length, in marked distinction to its ID counter- 
part Gi in ([8]). This independence characterizes the 2D 
trapped case as marginal [5| , in the sense that for any fi- 
nite trapping length (nonzero trapping frequency), single 
condensates persist as long as the interaction coupling is 
sufficiently small. Physically, this corresponds to the fact 
that the usual logarithmic divergence of (the integral of) 
phase fluctuations in two spatial dimensions 0] , is cut off 
by the trapping. 

The calculation then proceeds along the same lines 
as in the quasi- ID case and, again using ([7]), we have 
^NX - ^X'^. The continuum limit en- 
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ergy expression has a form analogous to Eq. 
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the relation G2UJ_l/{NAi) = A' 
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Note, however, the different power law dependence on A 
in the interaction terms. 



The minimum equation obtained by differentating (|16p 
with respect to A now is solvable analytically; it is a 
quadratic equation in A^, yielding 



(A2 



5G2 - 2 - Vl6-80G2 + (G2)2 (17) 



Due to the independence of G2 — ^J^j-KNasjlz on the 
trapping in the plane, condensate fragmentation in quasi- 
2D obtains at fixed transverse trapping solely because of 
an increase of the combination Na^. 



B. Three- mode approximation 



We now consider all three modes in p3)) as the inde- 
pendent field operator modes. *I'(r) — X]i=o ± ^«^*(^)- 
A significant simplification takes place in the Hamilto- 
nian due to the exact vanishing of the pair-exchange co- 
efficients involving two particles of the same mode ± and 
the single-particle ground state owing to the cylindrical 
symmetry of the excitation modes. The remaining pair- 
exchange process is between a pair of counter- and co- 
propagating azimuthal excitations and the ground state. 
The Hamiltonian then reads 

^ ^ 



J={±} 



(20) 



The interaction coefficient Bi = 2Trg J rdrdzip{ 
g/ii27rf /^BlUz), and 



B,^^, i3.M = |l,i,2,2^ (21) 



iVA2 ' 

^-lA^ 

£+ = £_ = 2eo- We evaluate the total energy in the 



The single-particle energies are Cg/wj, = 2^ + ^A"^ and 
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three-state Fock basis, with the following ansatz for the 
state vector 

|vI/)=^Vu'|A^-;-Z', /,/'), (22) 

Ll' 

that is with the probability amplitudes ipi^i' for I and /' 
particles being in the co- and counterpropagating modes, 



respectively. Note that the excited state in the two- 
mode approximation corresponds in the three-mode basis 
to setting 14-) = 4| E/ (1^^ - I, I. 0) + \N- I, 0, 1)). 



The ansatz ([22]) results in the following functional of 
the total energy 
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{mm = E iv^'-' 



eoN + (ei - 6o)a + I') + ^{N - I - l'){N - I - I' - I) + ^ [l{l - I) + - 1)] 



+ ^{N-l-l'){l + l') + ?±ll' 



B3 

+ —di,i'tpi i,tpi+i,i'+i + h.c. , 



(23) 



where the pair-exchange coefhcient here reads di^i' = 
[{N-l-l'-l){N-l-l'){l + l){l' + l)]^/^. We minimize 
with respect to the state amplitudes ipi^i' , to get 



B3 



(24) 



where q.z' = cqN + {e, ^ eo){l + 1') + ^{N - I - V){N - 
l-r-l) + ^ - 1) + /'(/' - 1)] + ^{N - I - l'){l + 

We argue that the peculiar matrix structure of ([M)) 
allows for an understanding of the [N -f \){N -\- 2)/2 
dimensional eigenvalue problem in {l^V) in the form of 
iV + 1 smaller problems indexed by fc = 0, . . . , [where 
k = — and only dependent on one running vari- 
able I. To see this, we rewrite the state amplitudes 
ifji^ii = '4'1,1+k = V'f ^nd also the matrix elements qj' , di^ii ; 
then, we note that only diagonal coupling terms occur in 
(HI, i-e. 
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(25) 



where the superscript k runs from to N, thereby re- 
stricting the corresponding I e {0, . . . , N/2~ [fc/2] }. The 
symmetry of the eigenvalue problem with respect to the 
exchange of I and I' justifies the reduction to I < I' that 
has been implicitly performed by choosing the ranges of 
k and /. Hence all states at a given fc 7^ are at least 
twofold degenerate. The above reformulation of ([M]) elu- 
cidates that different values of k belong to uncoupled 
equations and thus allows us to study the lower dimen- 
sional problems (j25p independently. In matrix language, 
this observation can be understood as a particular order- 
ing of the {I, I') indices, such that the eigenvalue problem 
decomposes into a block diagonal matrix, with each of 
the TV -|- 1 blocks corresponding to a symmetric tridiago- 
nal matrix for a given k by Eq. (|25|) and with block sizes 



N/2 + 1 — [fc/2] , where the square brackets indicate the 
next larger integer. 

In the large TV limit, we again apply the continuum 
limit to respecting the alternating signs for 

even and odd values of I in the , which occur due 
to the repulsive pair-exchange coupling ^3 > 0, cf. the 
corresponding distribution in the quasi-lD case, Eq. ([5|). 
Taking the continuum limit yields 

-^A|^^-| + imol}2(^_g)2|^,fe| ^ _ (26) 



with inverse mass coefficient of the harmonic oscillator 
analogy 1/mo = VW~fc)^(]vT3fc) , the effective 
frequency 



V2N^ 



'k^ + N^ 



8fc3 

N + 3k 



y/{N - ky^{N + 3k) 



(27) 



and the shift of the |'!/''^(0I distribution from / = N/A, 

G2W_L 



& = 



1 



Afk + N) + kJ ^^LJl 

4^ ^ V TV + 3fc 



A2 



(28) 



Finally, the I independent (but fc dependent) energy 
shift reads 



En = 



4A2 



^ {15N^ - 2kN - 5fc2) 



+ (1 + A^) (3iV + fc) 



1 

mo 



mon-'e\{29) 



In these expressions, we have substituted the matrix ele- 
ments Bi with their values obtained from the variational 
ansatz in (|2ip . and the corresponding single-particle en- 
ergies e,;. 

The full ground state energy in this limit reads E = 
51/2 -|- Eq and its variational minimum at Asm can be 
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FIG. 4. (Color online) Energy minima of the three-mode 
ansatz for the quasi-2D trapping geometry in terms of the 
variational parameter A, and three-mode degree of fragmen- 
tation 5 (inset) from Eq. (|32|l , both as a function of the in- 
teraction measure G2. Conventions are the same as those of 
Figs.Handm The number of particles here is iV = 20000. 



solved for analytically. We do not display the very 
lengthy resulting expression here. The Taylor expansion 
around k/N = reads, cf. the expression for the two- 
mode case in ([T7|) . 

2 

^4 - A4 , ^-2 (448 _ 151^2 



A'^ - A* ■ 



+0 



G2 

12G> 

3 



56G>)(- 



(30) 



with G> = ^16 -8OG2 + (6*2)2. 

Exactly as in the quasi-2D two mode case, at the tran- 
sition point where G> becomes real, given by (62)^, in 
P9I) . the continuum limit becomes valid, and matches 
the newly appearing minimum of the numerical solution 
of the eigenvalue problem (|25|) . We then precisely re- 
cover the two-mode energy in the limit iV — > 00 for 
fc = 0. This is a nontrivial result, because the corre- 
sponding states in the Fock basis do not coincide. Exam- 
ining the resulting energy minimum at Asm, we observe 
a monotonous increase with the eigenvalue index k and 
can hence justify focussing only on minimal energy solu- 
tions for states around k. Below, we discuss the effects of 
quasi degenerate states for small k whose occurrence can 
be understood in the continuum limit by showing that 

dkEmin{k)\k=0 = 0. 

For finite particle numbers, comparing the ground- 
state energies for the three- and two-mode ansatze, we 
find that the three-mode model yields a slightly lower 
energy at fc = 0, of order O.OOlwj^ per particle for par- 
ticle numbers of order N 10^, further decreasing with 
increasing N and asymptotically approaching zero due to 
the identity of the energies in the continuum limit. 

The continuum limit, as well as numerical experiments 
for finite N thus justify taking the minimal energy so- 
lution of the k ~ problem as the global ground state. 



which amounts to taking the state with equally populated 
counterclockwise and clockwise circulating single-particle 
modes. However, we note that because of the similarity 
in structure for the eigenvalue problems of neighbouring 
values of fc, the difference of the minimal energy eigenval- 
ues per particle corresponding to different k is of order 
uj_i_/N. For large particle numbers, we thus expect quasi- 
degcneracy of the eigenvalues located around fc = 0, be- 
cause they then become arbitrarily closely spaced. 

To study the consequences on fragmentation, we con- 
sider the one-particle density matrix whose off-diagonal 
elements are given by 

1,1' 

Y,ViN~i~m + i)^i,^i+iM, (31) 



l.V 



po- = J2 ViN~i~nii' + m,i,i^i,i'+i, 
1,1' 

and introduce a generalized fragmentation measure for 
three modes 



;?=l-|Ao-A2|/iV, 



(32) 



where Xi are the eigenvalues of the one-particle density 
matrix ordered by size with Ao being the largest eigen- 
value and A2 the smallest. 

Supposing that the system only occupies a single state 
A: = 0, it is apparent from the /c-diagonal form of (|25p 
that all off-diagonal elements of the density matrix in 
([31]) vanish. This is equivalent to a loss of coherence 
between states carrying different angular momenta. In 
the continuum limit, the fragmentation can be written in 



terms of the shift as = 1 



where the shift is 



^ N/3 

defined by (ajoo) = N/2 + 26 with -N/A < e < N/A 
(at k = 0). The fragmentation is then obtained to be, 
cf.Eq.dUD, 



5(G2) = 1 - 



I-G2/4 



G2 









VG2) 







(33) 



and has its asymptotic maximum value 33 % at large 
G2 3> (62)^, see also the inset of Fig.gl 

In the quasi-degenerate case of large N we expect many 
neighbouring low lying k states to be occupied, lead- 
ing for any small perturbation to a mixing of different k 
states away from A: = 0, and to nonvanishing off-diagonals 
(Plj) . Hereby coherence is established, and a nonfrag- 
mented ground state in terms of the three-mode single- 
particle density matrix and the corresponding fragmen- 
tation measure is obtained. On the other hand, super- 
position of many modes of different k leads to a (coher- 
ently superposed) excited state of zero angular momen- 
tum [taking into account both degenerate sectors corre- 
sponding to / < Z' and I > I'], and we are then led back 
to two-mode fragmentation between ground state and the 
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coherently superposed excited state. This is due to the 
fact shown in section IIII Al namely that the two-mode 
fragmented state is always lower in energy than a single 
condensate above a critical G2 given by ([T5)) . 

IV. DISCUSSION AND CONCLUSIONS 

We considered harmonically trapped low-dimensional 
gases and found by a variational analysis that they split 
from a single macroscopically occupied field operator 
mode in the weakly confining direction (the condensate) 
into two (quasi-lD trapping) respectively three (quasi-2D 
trapping) such macroscopically occupied modes with no 
remaining phase coherence between them, upon increas- 
ing a dimensionless measure of interaction strength be- 
yond a critical value. We have furthermore demonstrated 
that due to the symmetry of the matrix equations for 
the state vector amplitudes, three-mode fragmentation is 
highly susceptible to decay into a two-mode fragmented 
state in the quasi-2D isotropic trapping geometry. 

The results obtained represent, then, to the best of our 
knowledge the first example of ground-state fragmented 
scalar condensates in a single trap, and implement the 
result of in a physically realistic situation. They 
imply that the quasi-lD and quasi-2D condensates de- 
cay into fragmented condensates of a few macroscopi- 
cally occupied modes for a sufficiently large interaction 
energy, with the quasi-lD condensates, as could be ex- 
pected, being significantly more fragile. By way of a nu- 
merical example for the quasi-lD situation, trapping fre- 
quencies of iOzj^T^ = 3.5 Hz and lox^J^-k = 360 Hz (used 
in the magnetic trapping setup of [6|) correspond to a 
ratio lzll\. — 10. For ^^Na, we have /_l — 1 /im. For 
Os ~ 2.8 nm, the number of atoms driving G\ above its 
critical value is only of order ~ 500. This raises the 
interesting question what was actually observed in the 
experiment [6| , where the quoted particle numbers are of 
order N ^ 10^ for the quasi-lD setup. While the thresh- 
old to the first onset of fragmentation seems surprisingly 



low, one should note that the energy difference separating 
fragmented and nonfragmented condensate states is very 
small close to the threshold, Eq. ([TU)) . the energy barri- 
ers between the two states being even smaller, cf. Fig.[T] 
In order to observe the (zero temperature) transition 
point itself, extremely low temperatures would therefore 
be required. Fragmented condensates become energeti- 
cally clearly preferred over single condensates at experi- 
mentally accessible thermal fluctuation levels only at sig- 
nificantly higher values of the dimensionless interaction 
measures, and therefore at much higher particle numbers 
for given trapping frequencies. Another possibility to in- 
crease the interaction measures and to obtain fragmented 
condensate states is to increase the coupling constant it- 
self li]. 

We have shown that fragmented condensate states in 
low-dimensional systems can be obtained well before the 
thermodynamic limit of infinite extension. While the de- 
gree of fragmentation depends on the (number of varia- 
tional) modes chosen for a particular trapping setup as 
well as the applied fragmentation measure, the occur- 
rence of few-mode condensate fragmentation for inter- 
acting bosonic gases should be a rather generic feature 
of their many-body physics. Few-mode fragmentation is 
owed to the confined nature of the system and a suf- 
ficiently strong and positive pair-exchange coupling be- 
tween the bosonic modes. 
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